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ABSTRACT 



We offer a brief description of tlie bulk-motion Comptonization (BMC) model 
for accretion onto black holes, illustrated by its application to observational data 
for LMC X-1, and Nova Muscae 1991. We then extract some physical parameters 
of these systems from observables within the context of the BMC model, drawing 
from results on GRO J 1655-40, for which we presented extensive analysis 
previously. We derive estimates of the mass, (16 ± 1)Mq x [0.5/cos(i)]^'^ and 
mass accretion rate in the disk in Eddington units (m^; ^ 2) for LMC X-1, and 
(24 ± 1)Mq X [0.5/cos(i)]^/2^5.5 and (m^ ~ 3 ) for Nova Muscae 1991 [where 
cos(i) and (is. 5 are the inclination angle cosine and distance in 5.5 kpc units 
respectively]. Differences between these estimates and previous estimates based 
on dynamical studies are discussed. It is further shown that the disk inner 
radius increases with the high-to-low state transition in Nova Muscae 1991. 
Specifically, our analysis suggests that the inner-disk radius increases to r in — 17 
Scwarzschild radii as the transition to the low-hard state occurs. 

Subject headings: accretion — black hole physics radiation mechanisms: 

nonthermal — relativity — stars: individual (LMC X-1, Nova Muscae 1991) 
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1. Introduction 

It has long been known that accreting stellar-mass black holes in Galactic binaries 
exhibit a "bi- modal" spectral behavior - namely the so called high-soft and low-hard 
spectral states; see Liang (1998) for a recent review. High and low in this context allude to 
the relative 2-10 keV luminosities in a given system, which in turn responds to the rate at 
which mass accretion is taking place. An important point is that the high-soft state spectra 
are apparently a unique black hole signature, whereas the low-hard state spectral form has 
been seen in neutron star binaries under certain conditions - typically in low luminosity 
states of bursters (Barret & Grindlay 1995). It thus seems reasonable to speculate that 
this unique spectral signature is directly tied to the black hole event horizon. This is the 
primary motivation for the Bulk Motion Comptonization Model (BMC) introduced in 
several previous papers, and recently applied with striking success to a substantial body of 
observational data. We will not describe the model in detail here, but refer the interested 
reader to the literature and briefly state its essential components: Compton scattering of 
thermal X-rays from bulk-motion infall in close proximity to the black-hole event horizon, 
and a hot Compton cloud obscuring the inner disk during periods of low-mass-accretion rate 
[Titarchuk, Mastichiadis & Kylafis 1997, Titarchuk & Zannias 1998, Shrader & Titarchuk 
1998, Laurent & Titarchuk 1999, Borozdin et al. 1999 (hereafter TMK97, TZ98, ShT98, 
LT99, BOR99 respectively)]. 

Given the success of the BMC model in reproducing the high-soft state continuum, 
we have begun to investigate potential predictive power. From the spectral shape and 
normalization, one can calculate an effective disk radius and the mass-to-distance ratio 
m/d. Other quantities such as a black hole mass and mass-accretion rate can be determined 
if the distance or mass are known independently. This is dependent upon an additional 
factor, Thi the so called "hardening factor" which represents the ratio of color-to-effective 
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temperature. If one can identify several "calibrators", i.e. sources for which distance and 
mass are well determined, the model can be applied to a derive Th over the available 
dynamic range in luminosity. Three such sources are GRO J1655-40, for which we have 
presented extensive analysis in a previous paper, LMC X-1 for which the distance is well 
constrained, and Nova Muscae 1991 (herein NM91), for which the binary parameters are 
reasonably constrained. In the latter two cases, we have assumed the hardening factor, 
Th — 2.6, derived from our analysis of GRO J1655-40. The temperature-flux curves we 
present provide a self-consistency test of this assumption. 

In this paper, we apply the BMC model to LMC X-1 and NM91 then, drawing from 
previous analysis of GRO J1655-40, infer some physical parameters of these binary systems 
based on the model and on methodology detailed in BOR99 and section 3. In section 4 we 
offer interpretation of our results including qualitative inferences regarding the evolution of 
the emission-region size over the course of the NM91 decline. 

2. PHYSICAL PARAMETER ESTIMATION 

If the physics underlying the emergent spectrum is understood, one can, with minimal 
model dependency, use observables to determine parameters of the system such as the 
mass-to-distance ratio, the mass accretion rate, and the disk effective radius. From the 
inferred color temperature and absolute normalization, the effective area of the emission 
region can be obtained in terms of the distance and Th- The observable surface area can be 
related to the black hole mass using, for example, the Shakura-Sunyaev (1973), herein SS73, 
prescription modified to treat electron scattering and free-free absorption and emission 
[Titarchuk 1994 (hereafter Ti94), see also Shimura & Takahara 1995]. The emergent 
spectrum is then represented by the integral of a diluted blackbody function characterized 
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by Tf^'^. The flux density g^, measured in units of erg sec~^ ster~^cm~^, is given by: 

9AE) = C^-2 ) J" ' / rdr (1) 

Jr,„ exp[E/Tmaxj{r)\ - 1 

wliere r = R/Rg is in Scliwarzcliild and m = M/Mq in solar units. The color temperature 
distribution T^ax ■ f{r) which has the radial dependence, 

.. . ^ [1 - (r.Jr)VT/^(1.36r.03/4 

jy) ^3/4(1/7)1/4 ' y^) 

is normalized to T^ax at r = (7/6)^rj„ = 1.36rj„. The normalization factor of the equation 

(2) measured in units phot sec"^ ster^-*^ cm^ is Cn = 0.9l{m/ d)"^ cos iTf^'^ where i is the 

inclination angle (between the line of sight and the disk normal), and d is in kpc. We then 

apply the mean value theorem to replace the radial temperature distribution with T^oi 

giving 

, ^3(1.6- 10^9) 

exp[E/TcoiJ - 1 

where A^ = CN^'^ff- In general, re// is a function of the photon energy E, however if the 
energy coverage is such that the ratio E/kTcoi ^ 2, r^ff is nearly constant and the single 
temperature blackbody is a reasonable approximation to the disk spectrum. The value of 
Teff then depends only on T^oi (or Tmax)- BOR99 derived Th using inferred parameters and 
an averaging method over a large flux range [ (0.3 — 2.2) • 10^ ergs s~^cm^] obtaining that 
Th — 2.6. We can then calculate the mass-to-distance ratio 

(dV 0.91T,-Vg,, . 

— = ^cos? (4) 

\mj An 

and if the distance or mass are known, the mass accretion rate 

. _^ (T^..)1(7/6)V.„]3 
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3. VIABILITY OF THE METHOD: APPLICATION TO OBSERVATIONAL 

DATA 

To test the viability of our method, we have collected data for LMC X-1, for which 
the distance is well constrained by virtue of its location in the LMC - herein we adopt 
the LMC distance of 51.2 ± 3.1 kpc of Panagia et al (1991) as the distance to LMC X-1 
- and NM91 for which the binary parameters are reasonably well constrained and the 
distance is known to a factor of ~ 2. From earlier work on GRO J1655-40 (BOR99), we 
determined that T/j ~ 2.6, and demonstrated that the data are strikingly consistent with the 
Temperature-Flux relationship predicted by the SS73/Ti94 model. We use this information 
here. 

LMC X-1, has been observed extensively with various instruments, for example, 
HEAO-1 (White & Marshall, 1984), BBXRT (Schlegel et al. 1994), and RXTE (Schmidtke 
et al. 1999). Although luminous, given its ~ 50kpc distance, its intensity is much lower than 
that of typical Galactic XRBs. The impact on our study is a lack of reliable measurement 
above ~ 20 keV. We found that typical 2-10 ksec exposures available from the RXTE data 
archives did not include reasonable HEXTE (20-200 keV) measurements. We also examined 
the publicly available CGRO/OSSE data. The source appears to be marginally detected in 
a 200 ksec exposure, a simple power law fit yielding a spectral index consistent with typical 
PC A values. However, it is poorly constrained {a = 1.71 ± 1.2), and since the data were not 
contemporaneously obtained, we cannot reliably apply our model. Usually, however, the 
proportional counter data alone were sufficient to constrain our model parameters. NM91 
was a bright X-ray transient (also known as GS 1124-683 = GRS 1121-68) discovered in 
January 1991 independently by Ginga/ASM and GRANAT/Watch. It reached a peak 
intensity of about 6 times the Crab over a 2-10 keV band, and its light curve exhibited 
the canonical X-ray nova fast-rise exponential-decay profile (e.g. Chen, Shrader & Livio 
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1997). Subsequent to its discovery, a number of follow-up observations with the Ginga LAC 
detectors were made, covering portions of the outburst peak and subsequent decay. Both 
the high-soft and low-hard states were seen (Ebisawa et al. 1994). 

We then followed procedures detailed in BOR99 to analyze 7 spectra for NM91 and 
about 30 spectra for LMC-Xl. Some representative results, for LMC X-1 (averaged over 
the sample) are T=0.79 ± 0.22, a = 1.63 ± 0.24 and / = 0.17 ± 0.23 (also see Figures 1 
and 2.). Cases of particular interest for NM91, are the March 1991 data (/ ~ 0) when the 
hard-power- law component is not detected, and September when the source has entered 
the low-hard state (T,q;, /) = (0.29, 0.89, 0.37). A typical peak-outburst high-state spectrum 
(e.g. 20 January 1991) has {T,a, /) = (0.81, 1.36,0.15). While we were unable to obtain a 
satisfactory fit to the March 1991 data, we found that we could model the September 1991 
data not only in terms of a thermal Comptonization spectrum (or simple hard power law), 
but as a BMC spectrum with a relatively large /—value - i.e. the scattered component 
had become dominant. The disk at this stage is nearly obscured and our spectral index 
parameter decreases to values below the asymptotic limit of the typical high-state case 
(TZ98). 

The apparent broad absorption feature, first noted by Ebisawa et al. (1994), and 
modeled in detail by Zycki, Done & Smith (1998) is evident in our deconvolved spectrum 
of 12 January, 1991. However, the depth that we obtain is only about half that of Ebisawa 
et al. (1994), and it is substantially less pronounced in the 20 January 1991 spectrum. We 
do not find strong evidence for such a feature in subsequent spectra, however, statistics 
become limiting as the source faded. In order to treat the entire evolutionary sequence in a 
consistent manner, we have opted not to incorporate additional features into our analysis. 

Following the prescription outlined in the previous section, we have constructed a plot 
of the empirical fiux-temperature relationship for LMC X-1 and NM91 (Figure 3). The 



solid curves represent the corresponding relationship in the context of our disk model for 
the indicated value of rj„, while the dashed-line curves represent the 20% deviation corridor 
about the model for rj„ = 3. This utilizes the hardening factor obtained by BOR99. The 
theoretical curves can be reproduced analytically using the following formulae: 

F = rlffT^oi / , -^^dx. (6) 

We find that r^ff can be represented by a fourth-order polynomial of the form anT^oi 
with the coefficients ao = 2.674, 4.452, 6.737, 14.33, ai = 14.09, 24.57, 43.61, 96.07, 
02 = -13.78, - 30.51, - 66.11, - 170.1, a^ = 14.75, 40.76, 94.68, 282., 
04 = -5.42, - 19.77, - 49.5, - 176.1 for n^ = 3, 5, 8, 17 respectively. We 
then applied a normalization factor of 10.1 to the theoretical values to match the data. 

The data in Figure 3 were obtained by integrating our model flux from 4 — 20 keV. 
A scale factor of 200 was applied to LMC X-1 to present both sources in the same plot 
(200 is approximately the square of the ratio of the mass-to-distance values we derive for 
NM91and LMC X-1). Most of the points are associated with the inner disk edge near 3Rs. 
However, as the flux decreases, and NM91 enters the low-hard state, the inner edge recedes 
outwards to 17 Schwarzchild radii. Esin et al. (1997; hereafter EMCN) argue that the disk 
is located near the last stable orbit at rj„ = 3 in the high state, consistent with what we 
infer. However, our conclusions are incompatible with the EMCN scenario of an accretion 
flow extending out to rj„ ~ 10^ during decline, requiring instead a transition from the high 
to low states and optically thin to optically thick structure at rin — 17. Zycki et al. (1997) 
have independently reached the same conclusion, that the retreat of the inner disk during 
the high to the low state transition is much slower than predicted by EMCN and optically 
thick material within ~ 10 — 12 i?^ is generally present. The path we have followed in 
reaching these conclusions however, is distinct from Zycki et al (1997), primarily in that 
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those authors applied an additional model component - an absorption-edge feature - from 
which information regarding the disk was extracted. We have applied a self-consistently 
calculated model without additive components. Using extensive hydrodynamical and 
radiative transfer calculations, Chakrabarti & Titarchuk (1995) first predicted, that the 
apparent position of the inner edge of the accretion disk would be seen as moving to the 
position of the shock location (i.e. between 8 and 15 Rg). Our conclusions here would seem 
to support that prediction. 

4. DISCUSSION AND CONCLUSIONS 

We have calculated mass-to-distance ratios following the method outlined in the 
previous section. For NM91, Orosz et al. (1996) have presented a comprehensive dynamical 
study of the quiescent binary system, leading to a maximum probability contour in 
m — sin(i) space with m ranging from about 5 to 7.5 Mq for inclination i ~ 55° — 65°. The 
most uncertain parameter however is its distance, with published estimates ranging from 
1.4 kpc (Delia Valle et al 1991) to 8 kpc (Cheng et al. 1993). Orosz et al. (1996) who argue 
that d = 5.5 ±1 kpc based on their characterization of the secondary star, and brightness 
estimates for the quiescent disk, however, their uncertainty may be understated. We 
calculate the BH mass, m = [24 ± 1] x [1/2/ cos{i)Y^'^drj,5(Th/2.6Y and mass accretion rate 
in the disk in Eddington units, rrid ~ 2.5 [using formula (5), mass value 24 and Tmax ~ 0.95 
keV] where ^5.5 = d/5.5 kpc. Our mass estimate can be reconciled with Orosz et al. (1996) 
only if (^55 < 0.31, which would make the outburst luminosity ^3-20 — 1-5 x lO^^ergs s~^. 
We note, that for our fiux measurement, if ^5.5 = 1 and the mass is near the lower end of 
the range cited by Orosz et al (1996), the luminosity is very nearly at the Eddington limit, 
whereas in our case it is substantially sub-Eddington. However, since stable disk accretion 
is occurring, super-Eddington levels cannot persist and a mass as low as ~ 5Mq seems 
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unlikely. 

For LMC X-1, we can calculate the mass directly. We find m = 
[16 ± 1] X [(0.5/cos(i)]^/^(T/i/2.6)^ and m^ ^ 1.7. This large of a mass is consistent 
with the fact that LMC X-1 is relatively luminous compared to typical Galactic X-ray 
binaries. It would seem, however, to be at odds with published values obtained from 
radial velocity measurements of the candidate optical companion "star 32", of 4 — IOMq 
(Hutchings et al 1987). To reconcile our result with the upper allowed-range of their 
analysis would require a low-inclination system, i < 30°, however, this could render 
the ellipsoidal-variability amplitude below the detectability threshold. It must be noted 
however, that the association of LMC X-1 with star 32 remains speculative, and the only 
significant constraint on the inclination is an apparent lack of eclipsing. We can make an 
additional argument against the lower end of the Hutchings et al (1987) mass range: for 
a given disk inner radius and temperature distribution, tti scales in proportion with m. 
Since m/d, and in this case m is an invariant of our model, for m < 10, m < 1, which is 
inconsistent with a high-soft state source (LT99, BOR99). 

The temperature flux relation in NM91 seems to populate two distinct tracks. This 
is an indication that the emission area, and thus the inner-disk radius, is increasing. This 
is consistent with the previous conclusions of Zycki et al. (1998) who find rj„ ~ 13 - as 
illustrated in Figure 3 an inner disk radius of rj„ ~ 17 is consistent with our analysis. 

An issue central to stellar evolution theory is the distribution of black hole masses in 
the Galaxy and contact binaries may provide our only insight. Dynamical analyses based on 
optical studies are often limited due to dust obscuration or the brightness dominance of the 
disk. Our method offers, in principle, an independent means of parameter determination. 

ACKNOWLEDGMENTS 



-11- 



This work made use of the High-Energy Astrophysics Science Archive Research Center, 
the Compton Gamma Ray Observatory and Rossi X-Ray Timing Explorer science support 
facihties at the NASA Goddard Space Fhght Center. We also wish to acknowledge Ken 
Ebisawa for providing the Ginga data and for useful discussions. 



- 12- 

REFERENCES 

Barret, D., & Grindlay, J.E. 1995, ApJ, 440, 841. 

Borozdin, K., Revnivtsev, M., Trudolyubov, S., Shrader, C.R., & Titarchuk, L.G. 1999, 
ApJ, 517, tbd. (BOR99) 

Chakrabarti S.K. & Titarchuk, L. G. 1995, ApJ, 455, 623 

Chen, W., Shrader, C.R., & Livio, M. 1997, ApJ, 491, 312 

Cheng, F.H., Home, K., Panagia, N., Shrader, C.R., Gilmozi, R., Paresce, F., & Lund, N. 
1992, ApJ, 397, 664. 

Hutchings, J.B., Crampton, D., Cowley, A. P., Bianchi, L., & Thompson, LB. 1987, A J, 94, 
340. 

Delia Valle, M., Jarvis, B.J., & West, R.M. 1991, Nature, 353, 50. 

Ebisawa, K., et al. 1994, PASJ, 46, 375 

Esin, A.A., McClintock, J.E., & Narayan, R. 1997, ApJ, 489, 865 (EMCN) 

Laurent, Ph., & Titarchuk, L.G. 1999, ApJ, 511, 289 (LT98) 

Liang, E.P. 1998, Phys. Rep., 302, 66. 

Orosz, J.A., Bailyn, CD., McClintock, J.E., & Remillard, R. 1996, ApJ, 468, 330. 

Panagia, N., Gilmozzi, R., Macchetto, F., Adorf, H.-M., & Kirshner, R. P. 1991, ApJ, 1991, 
380, L23 

Schmidtke, P.C, Ponder, A.L., & Cowley, A.P. 1999, ApJ (astro/ph-9812358) 

Schlegel, E.M., et al. 1994, ApJ, 422, 243. 



-13- 

Shakura, N.I., & Sunyaev, R.A. 1973, A&A, 24, 337 (SS73) 

Shimura, T., & Takahara, F. 1995, ApJ, 445, 780 

Shrader, C.R., & Titarchuk,L.G. 1998, ApJ, 499, L31 (ShT98) 

Titarchuk, L.G. 1994, 429, 340 (Ti94) 

Titarchuk, L., & Zannias. T. 1998, ApJ, 493, 863 (TZ98) 

Titarchuk, L. G., Mastichiadis, A., Kylafis, N. D. 1997, ApJ, 487, 834 (TMK97) 

White, N.E., & Marshall, F.E., 1984, ApJ, 281, 354 

Zycki, P., Done,C., & Smith, D. 1998, ApJ, 496, L25. 



This manuscript was prepared with the AAS h^T[T;X macros v4.0. 



-14- 



Fig. 1. — Typical BMC model- fit to the LMC-Xl high-energy continuum, with the inferred 
parameters in the upper right. Data are from RXTE/PCA observations on 16 April, 1997. 
The residuals are normalized to the model values. We found the source intensity varied by 
about a factor of 3 and Tqol ranged from about 0.7 to 0.9 over our data sample. 



Fig. 2.— Typical BMC model-fit to the NM91 high-energy continuum. Model-fit to NM91 
high-energy continuum which yielded the indicated parameters. These data resulted from 
a 20 January, 1991 observation with the Ginga/LAC when about 10 days after the initial 
outburst 



Fig. 3. — Temperature-fiux relation for NM91 and LMC-Xl. The solid curves represent 
the model predictions for the indicated inner-disk radii, calculated according to the method 
presented in the text. The dashed curves represent a 20% deviation about the model. The 
points which lie on the lower curve, corresponding to an inner radius of 17, are from the 
low-state data for NM91. 
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